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ABSTRACT 

We use the 2dF Galaxy Group Catalogue constructed by Merchan & Zandivarez to 
study the large scale structure of the Universe traced by galaxy groups. We concentrate 
on the computation of the power spectrum and the two point correlation function. The 
resulting group power spectrum shows a similar shape to the galaxy power spectrum 
obtained from the 2dF Galaxy Redshift Survey by Percival et al., but with a higher 
amplitude quantified by a relative bias in redshift space of 6s (fc) 1.5 on the range of 
scales analysed in this work, 0.025 < /c/ZiMpc""'^ < 0.45. The group two point corre- 
lation function for the total sample is well described by a power law with correlation 
length So = 8.9 ± 0.3 /i-^Mpc and slope 7 = -1.6 ± 0.1 on scales s < 20 /i"^Mpc. 
In order to study the dependence of the clustering properties on group mass we split 
the catalogue in four subsamples defined by different ranges of group virial masses. 
Our results are consistent with a 40% increase of the correlation length sq when 
the minimum mass of the sample increases from > 5 x 10^^ h~^MQ to 

M.vir > 1 X 10^** h~^MQ. These computations allow a fair estimate of the relation de- 
scribed by the correlation length sq and the mean intergroup separation dc for galaxy 
systems of low mass. Our results show that an empirical scaling law sq — 4.7 dj^'^^ 
provides a very good fit to the results from this work, as well as to previous results 
obtained for groups and clusters of galaxies. The same law describes the predictions for 
dark matter haloes in N-body simulations of ACDM models. We also extend our study 
to the redshift space distortions of galaxy groups, where we find that the anisotropics 
in the clustering pattern of the 2dF group catalogue are consistent with gravitational 
instability, with a flattening of the redshift-space correlation function contours in the 
direction of the line of sight. The group pairwise velocities found from this analysis 
for a sample of groups with masses Myir > 5 x 10^^ H^^Mq, are consistent with 
^y;2^i/2 _ (280^^5o)^'^/^' agreement with ACDM cosmological simulations. The 
bias factor for the 2df groups of moderate masses A4yir > 2 x 10^^ h~^MQ is consis- 
tent with the values predicted by the combination of a CDM model and the ellipsoidal 
collapse model for the formation of structures. 

Key words: galaxies: clusters: general - cosmology: large-scale structure of Universe. 



1 INTRODUCTION 

The position of groups of galaxies in the structure hierarchy 
makes them very interesting objects, as the density fluctua- 
tions sampled by groups lay between those traced by galax- 
ies and clusters of galaxies. This is what makes the study 
of galaxy groups a key area of research in cosmology and 
galaxy formation. Much efi'ort has been devoted to the study 
of groups in order to understand the large scale structure of 
the Universe (e.g., Jing & Zhang 1988, Maia & da Costa 
1990, Ramella, Geller & Huchra 1990, Trasarti-Battistoni, 



Invernizzi & Bonometto 1997). The main tool adopted in 
most of these works for the study of the spatial distribution 
of groups is the redshift space two point correlation function 
^(s). These studies show, as expected from the hierarchical 
model, that the amplitude of the group ^(s) falls roughly 
between the corresponding amplitudes of correlations ob- 
tained from galaxies and cluster of galaxies. Using a suit- 
able statistical sample of groups, Merchan, Maia & Lambas 
(2000) analysed the dependence of the two point correlation 
function on mass, by splitting the group sample in different 
ranges of virial mass, and extending this analysis to lower 
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mass systems (See also Bahcall & West, 1992, Croft et al. 
1997, Abadi, Lambas & Muriel, 1998, Borgani et al. 1999, 
Collins et al. 2000). Their results show that the amplitude of 
^(s) tends to increase significantly with the mass of the sam- 
ple. They measure a 50% increase in correlation length for 
a sample with mean inter-group separation dc ^ 28/i~^Mpc 
with respect to a sample with dc = 9?i~^Mpc. 

On the other hand, the measurement of the power spec- 
trum of density fluctuations, a widely used statistical tool 
for galaxies and clusters of galaxies, has seldom been ap- 
plied to the study of galaxy groups. This could be due to 
the fact that up to now, a statistically reliable sample of 
galaxy groups, from which a robust determination of the 
power spectrum can be made, has not been available. 

Other important constraints on cosmological models 
can also be obtained from the study of the anisotropies in the 
redshift-space correlation function. The three-dimensional 
clustering information available in rcdshift surveys is sub- 
ject to the influence of peculiar velocities on the distance 
measurements. This induces an anisotropy into the cluster- 
ing pattern, which can be used to determine cosmological 
parameters due to the fact that the amplitude of peculiar 
velocities is related to the matter density, and the linear 
bias parameter, h. This study is carried out by measuring 
the apparent distortion of the clustering pattern in the two 
point correlation function in redshift space. On large scales, 
where the linear perturbation theory approximation for the 
growth of density fluctuation is valid, flows of mass from 
low to high density regions generate an artificial enhance- 
ment of the density contrast so as to produce a compression 
of the ^ contours along the linc-of-sight direction (Kaiser 
1987). This phenomenon allows estimates of the parameter 
P ~ QP'^ /h to be made. On small scales, the main effect is 
to reduce the amplitude of density fluctuations by means of 
random non-linear motions in virialized regions, producing 
the "Fingers of God" effect. In this regime, it is possible 
to obtain estimates of the one-dimensional pairwise rms ve- 
locity dispersion (Davis & Peebles 1983). Such an analysis 
was recently applied to groups of galaxies identified in the 
Updated Zwicky Catalogue by Padilla et al. (2001). They 
found a one-dimensional pairwise rms velocity dispersion 
for groups of 250 ± 110km s~^ and a noisy estimate of 
the /3-parameter, (3 < 1.5. These results suggest the need 
for larger samples of groups in order to obtain more reliable 
results from these type of objects. 

The 2dF Galaxy Group Catalogue (hereafter 2dFGGC), 
constructed by Merchan & Zandivarez (2002) is one of the 
largest group catalogues at present. These groups were iden- 
tified from the 2dF public lOOK data release using a modified 
Huchra & Geller (1982) group finding algorithm that takes 
into account the magnitude limit of the 2dF sample, red- 
shift completeness, and angular masks (Colless et al. 2001). 
This sample of groups has been used to study the global 
effects of group environment on star formation (Martmez et 
al. 2002a), the effect of local environment on galaxy member 
spectral types (Dommguez et al. 2002), a statistical analy- 
sis of luminosity functions (Martmez et al. 2002b) and a 
compactness analysis of groups (Zandivarez et al. 2002). 

In this work we use the 2dFGGC as a tracer of the 
large scale structure of the Universe and compute the power 
spectrum of galaxy group fluctuations and the two point 
correlation function. The sample is being splitted by taking 



all galaxies above different cuts in virial mass in order to 
study the dependence of the clustering strength on group 
mass. Finally, we study the redshift-space distortions of the 
group two point correlation function and obtain estimates 
of the group pairwise velocities and the group bias factors. 

The layout of this paper is as follows. In section 2 we 
describe the group catalogue, while in section 3 we present 
detailed information about the construction of the mock cat- 
alogues used in the clustering analysis. The methods for es- 
timating the power spectrum and the two point correlation 
function using the 2dF group sample are described in section 
4. Studies of the anisotropies present in the redshift-space 
correlation function arc analysed in section 5. Finally, we 
summarise our conclusions in section 6. 

2 THE 2DFGGC 

Merchan & Zandivarez (2002) identified galaxy groups in the 
2dF public lOOK data release of galaxies within the south- 
ern (SGP, -37°. 5 <&< -22°. 5, 21'' 40" < q < 3''30'") and 
northern (NGP, -7°.5 < 5 < 2°.5; 9''50'" < Q < 14''50™) 
strips of the catalogue (see Colless et al. 2001). The group 
finder algorithm used in the identification is an adaptation of 
the algorithm developed by Huchra & Geller (1982), modi- 
fied in order to take into account the incomplete sky coverage 
of the 100k release of 2dF galaxies. 

The 2dFGGC was constructed using density contrast 
of 5p/ p = 80 and a fiducial linking length velocity of 
Vb = 200km s~^ which maximise the group finding accu- 
racy (see section 4 of Merchan & Zandivarez 2002). The 
resulting group catalogue contains systems with at least 4 
members, mean radial velocities in the range 900km s~^ < 
V < 75000km s~^ and a total number of 2198 groups (see 
Martmez et al. 2002b). The limit adopted in the number 
of members in galaxy groups is necessary in order to avoid 
pseudo- groups. 

The virial group masses, Mvir, are estimated us- 
ing virial radii, Rv, and velocity dispersions, a, M.vir = 
a^Rv/G. The Rv are computed using the projected virial 
radius and the a using the virial radius along the lino of 
sight (see section 5 in Merchan & Zandivarez 2002). A ro- 
bust estimation of this component is obtained by applying 
the bi-weight estimator for groups with richness Ntot > 15 
and the gapper estimator for poorer groups (Beers, Flynn 
and Gebhardt 1990, Girardi et al. 1993, Girardi and Giuricin 
2000). These methods improve the velocity dispersion es- 
timation in terms of efficiency and stability when deal- 
ing with small groups. The group catalogue has an av- 
erage velocity dispersion of 265 km s~^, an average virial 
mass of 9.1 x 10^^ Mq and an average virial radius of 
1.15 ft-^Mpc. 

The analysis involving computation of comoving dis- 
tances are carried out adopting the cosmological model pa- 
rameters Q.Q = 0.3 and = 0.7. 



3 MOCK CATALOGUES 

Throughout this work we use a set of mock catalogues con- 
structed from several cosmological numerical simulations of 
flat, low density, cold dark matter universes. We perform 
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Figure 1. Distance distributions of 2dF 100k release galaxies 
and groups. The left panels show the normalised distance distri- 
butions of galaxies in the 2dF galaxy redshift survey (2dFGPi,S) 
splittcd in the northern (NGP, upper panel) and southern (SGP, 
lower panel) strip. The centre panels show the normalised dis- 
tance distributions for the group catalogue (2dFGGC), while the 
right panels show the radial density distributions for the same cat- 
alogue. The solid lines in centre and right panels show the best 
fitting functions obtained for the group distance distributions (see 
equation 6). 



these simulations using the gravity part of the Hydra N- 
body code developed by Couchman et. al (1995), with 128^ 
particles in a cubic comoving volume of 180 Mpc on 
a side starting at z=50. The adopted cosmological model 
corresponds to a universe with a present day matter den- 
sity Qm = 0.3, vacuum energy density Oa ~ 0.7, baryon 
density fit = 0.0, spectral slope n = 2, T = 0.21, Hub- 
ble constant Ho = 100/ikm Mpc"'^ with h = 0.7, and 
an amplitude of mass fluctuations of crs = 0.9. In order to 
reproduce the same radial distribution as in the 2dFGRS, 
we adopt a galaxy luminosity function fitted by a Schechter 
function with M;^ - Slogig/i = -19.66, a = -1.21, = 
1.68 X 10~^ft^Mpc~^ and a model for the average k-|-e cor- 



rections given by the formula 



k{z) + eiz) 



(1) 



1 + 20«3 

(Norberg et al. 2002). Consequently, we assign absolute mag- 
nitudes to the simulation particles in order to reproduce 

the previous luminosity function and compute the apparent 
magnitude for the particles lying within the 2dFGRS mask, 
using the equation 

m = k + e + 5 logio(di/;i"^Mpc) -h 25 -h (M;, - 5 logio 

Then we apply the magnitude limited cuts using the mag- 
nitude limited mask of the 2dFGRS. In order to identify 
groups in the mock catalogues, we adopt the same values 
Sp/p = 80 and Vb = 200 km which maximise the group 



Figure 2. Upper panel: The comparison between the power spec- 
trum of the groups in the 2dFGGC (circles) and the power spec- 
trum of the galaxies in the 2dFGRS 100k release (triangles). The 
error bars are computed measuring the dispersion over 10 mock 
catalogues constructed from N-body simulations with a ACDM 
cosmology. Lower panel; The redshift space relative bias bs{k) be- 
tween galaxies and groups in the 2dF survey. The error bars are 
computed using the usual formula of error propagation. 



finding accuracy on the 2dFGGC. It should be stated that 

the resulting power spectra obtained with the mock group 
catalogues following the previous recipes are very good pre- 
dictions of the observed one. 



4 MEASURING THE SPATIAL 
DISTRIBUTION 

4.1 The power spectrum 

When calculating the power spectrum of the 2dFGGC we 
choose the method based on the Fast Fourier Transform 
(FFT) developed by Feldman, Kaiser & Peacock (1994, here- 
after FKP) in the version described by Hoyle et al. (1999). 
We denote by Ng the number of groups in the catalogue, 
and use the vector positions Xc to specify their positions. 
The following step in the FKP method is to construct an 
unclustered catalogue with the same geometry and radial se- 
lection function than the group sample. In order to minimise 
the contribution of Poisson errors to our uncertainties, this 
random catalogue should have a large number of points A^r, 
with vector positions x^. We define a quantity with mean 
value zero as: 

J(k) = uj{xc) e* "= - a e*""- (3) 

c— 1 r — 1 

where the first term on the right hand side is the Fourier 
transform of the spatial distribution of groups, and the sec- 
ond term is the Fourier transform of the window function 
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of the survey. Assuming Gaussian density fluctuations, FKP 
derive a weighting function 



oj{r) 



l+n(r)P™(fc)' 



(4) 



that minimises the power spectrum variance, where n{r) 
is the mean radial density of the catalogue and Pm{k) is 
the power spectrum. Given that in order to compute these 
weights we need an estimate of the power spectrum, we pro- 
pose a constant value for P„(fc) as an initial guess and then 
test the dependence of the results on this value. However, 
the result of a series of computations of the power spectrum 
varying this initial guess indicate that the assumed value of 
Pw{k) is not critical. A more difficult issue is the Eissumed 
mean radial density for our catalogue. As has been stated 
in many works on the 2dFGRS, the 100k public release of 
the catalogue has a completeness mask that should be taken 
into account for a reliable n(r) estimation. lu the computa- 
tion of the power spectrum from the 2dF galaxies, Percival 
et al. (2001) dealt with this problem by adopting a mean 
radial density which depends on the position on the sky. In 
order to search for a possible variation of the radial distri- 
butions with the choice of 2dF strip of the catalogue, we 
show in Figure 1 the normalised distance distribution of the 
2dFGGC in the northern (middle upper panel) and south- 
ern (middle lower panel) strips. The left panels show the 
same distribution for the galaxies in the 2dFGRS. The ra- 
dial densities for the group catalogue are plotted in the right 
hand side panels also for both strips as in previous panels. 
From these histograms one can appreciate slight differences 
in the radial distribution of groups in different strips. Thus, 
in order to take into account this variation, we use different 
mean densities for the two strips in the group catalogue. 



ni{r,a,S) 
n2{r,a,5) 



NGP 
SGP 



(5) 



where the angular dependence (a, 5) is related to a correc- 
tion that accounts for the redshift completeness mask of 
the 2dFGRS. This double prescription is also adopted for 
the construction of the random catalogues in order to im- 
pose more realistic radial selection functions. We use smooth 
curves which are the best fitting functions that describe 
the histograms showed in Figure 1. These functions are de- 
scribed by 



N{r) 



Nrr 



(^)'h(^)' 

V ' max / L V * max / 



-{1 + x/y) 



(6) 



where the parameters axe estimated using a chi-squaxe max- 
imum likelihood method. The fitting functions are plotted 
with solid lines in central and right panel of Figure 1. The 
best fitting parameters are x = 1.02, y = 7.13, Nmax = 71.44 
and Tmax = 362.34 7i"^Mpc for the NGP strip, and x = 1.03, 
y = 7.98, Nrr,ax = 66.20 and Vmax = 413.50 /i"^Mpc for the 
SGP strip. 

In Equation 3 we use a = Sc/Sr, where Sc = 
^^^^ w^(a;c) and Sr — '^^1^ Lj^{xr) in order to recover the 
definition of P(k) given by equation (3.4.5) of FKP. Then, 
the power spectrum estimator is obtained by: 



P(k) 



V{\5{k)\^ -ail + a)Sr) 



(7) 



n — I — I — r 
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s„=(8.9±0.3) h-' Mpc 
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Figure 3. The two point correlation function for groups in the 
2dFGGC. The error bars are estimated measuring the disper- 
sion over 10 mock catalogues as previously applied in the power 
spectrum error computation. The solid line is the best power law 
approximation (eq 9, see legend). 



where V is the volume over which periodicity is assumed. 
Finally, assuming isotropy we compute the power spectrum 
estimator averaging over spherical shells k < |k| < k + dk 
in A;-space. We compute spectral densities at multiples of 
the fundamental mode in order to avoid oversampling of the 
spectrum that could result in spurious features. The com- 
putation of these quantities is carried out by embedding the 
distributions within a periodic cubic volume V = r^„^ which 
is larger than the observational sample, and is divided into N 
cells per side with the spatial distribution of points (groups 
or random) assigned to the grid by means of the nearest grid 
point weight assignment scheme. 

The fiducial values chosen for the parameters involved 
in the calculation of the power spectrum are the FFT grid 
dimension A'' = 256, the side of the box, rtox = 2rmin, 
where rmin = 1244 /i~^Mpc is the side of the minimum box 
that contains the total catalogue (i.e., Nyquist frecuency 
~ 0.32 h~^Mpc^), a constant value of the power spectrum 
for the weight function of P„(fe) = 20000 /i"^Mpc^, and 
a number of random points of Nr = 1 x 10®. The resulting 
power spectrum for the 2dFGGC is plotted with filled circles 
in the upper panel of Figure 2. The error bars of this esti- 
mate are computed using a set of 10 mock catalogues con- 
structed as described in previous section. We also estimate 
the power spectrum for the galaxies in the 2dF survey with 
fiducial parameters N = 512 and P„(fc) = 5000 /i"^Mpc^, 
obtaining very similar results to those obtained by Percival 
et al. (2001). Our estimate of the galaxy power spectrum 
is shown by filled triangles in the upper panel of Figure 
2, where the error bars are also estimated from mock cat- 
alogues. As expected (Padilla & Baugh 2002), the power 
spectrum of galaxy groups has a very similar shape to the 
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Figure 4. The power spectrum for different mass subsamples 
of the 2dFGGC. The sohd hne shows the power spectrum for 
groups with A4vir > 5 X lO^'^h^^ Mq; dotted line shows the 
power spectrum for groups with M^ir > 2 X Mq; the 

dashed line displays the power spectrum for the subsamplc with 
Mvir > 5 X lQ^^h~^ Mq, and the dashed-dotted line shows the 
corresponding power spectrum for the subsample in the range 
Mvir > 1 X IQI^/i-IMq. The filled circles show the power spec- 
trum for the whole sample. 



result for galaxies but with a higher amplitude. This differ- 
ence in amplitude can be measured computing the relative 
bias between the two power spectra defined by 



Pgrpjk) 
Pgalik)- 



(8) 



The bias function is plotted in the lower panel of Figure 2 
where the error bars are computed from error propagation. It 
can be seen that the redshift-space bias parameter is almost 
constant bs{k) ~ 1.5 for the full range of wave numbers, 
0.025 < A;//iMpc~^ < 0.45. It should be taken into account 
that this redshift-space bias function could differ from the 
real-space bias function. A possible relation between them 
will be discussed in section 5. 



4.2 The two point correlation function 

As a second method to measure the clustering properties of 
the galaxy groups in the 2dFGGC we estimate the redshift- 
space correlation function. This computation is carried out 
using the estimator proposed by Landy & Szalay (1993) 



DD - 2DR + RR 
RR ' 



(9) 



where DD, DR and RR are the suitably normalised num- 
ber of data-data, data-random and random-random pairs re- 
spectively, in each separation bin. To estimate the two point 
correlation function, we generate a catalogue of randomly 



Figure 5. The two point correlation function for different mass 
subsamples of the 2dFGGC. The filled circles show the two point 
correlation function computed using the Landy & Szalay (1993) 
estimator (see equation 9). The dashed line corresponds to the 
best power law fit to the jjoiiits (sec labels), and the solid lines 
show the two point correlation function obtained from the inte- 
gration of the corresponding power spectrum (see equation 11). 



placed points with the same angular and radial selection 
function as the real data, as used for the computation of the 
power spectrum. We also take into account the completeness 
of the survey as a function of the position on the sky. 

The resulting two point correlation function for groups 
in the 2dFGGC, plotted in Figure 3, shows a positive signal 
up to 20 h~^Mpc. We adopt a power law approximation for 
our estimate described by the formula 



C(s) 



(10) 



where so is the correlation length and 7 the logarithmic slope 
of the correlation function. From a Levenbergh-Marquardt 
method (Press et al. 1986) which takes into account errors 
and applies a minimum nonlinear least-squares, the best- 
fitting parameters obtained for our estimation are so = 8.9± 
0.3 /i~^Mpc and 7 = — 1.6±0.1. This result is consistent with 
previous estimates obtained from samples with comparable 
space densities, by Girardi, Boschin & da Costa (2000) (so — 
8±1 h~^Mpc; 7 = -1.9±0.7) and Merchan, Maia & Lambas 
(2000) (so = 9.0 ± 0.4 ft-^Mpc; 7 = -1.67 ± 0.09). 

To test the dependence of the clustering of the 2dFGGC 
on the space density of our group samples, we study a set 
of subsamples defined by different ranges of virial masses 
which translate into values of mean inter-group separations, 
dc- The adopted mass limits are 

Sample 1 — » Mvir > 5 x 10^^ H'^Mq 
Sample 2 — » Mvir > 2 x lo" h~'^ Mq 
Sample 3 — » Mvir > 5 x 10^^ HT'^Mq 
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obtained from a minimum nonlinear least-squares analysis 
are 



40 60 

dyh-'Mpc 

Figure 6. Correlation length so as a function of the mean in- 
tergroup separation dc- Filled circles correspond to the four sub- 
samples defined in the 2dFGGC for different ranges of virial mass. 
Open circles show the determinations for subsamples of groups in 
the GUZC (Merchan, Maia & Lambas 2000). The open squares 
show the So — dc relation of APM clusters, and crosses correspond 
to the prediction for a A CDM model, both determined by Croft 
et al. (1997). The dashed line represents the universal scaling law, 
so = 0.4 dc (Bahcall & West 1992). The solid Une is our fit to the 
data, so = 4.7 dc°-^^. 



Sample 4 



Mvir > 1 X 



10" h- 



The analysis of these subsamples is carried out as follows. 
We compute the power spectrum for each subsamplc. These 
determinations can be considered to be statistically reliable 
due to the large number of groups in the 2dFGGC. Then, 
we make an estimate of ^(s) for each subsamplc using the 
relation between the power spectrum and the correlation 
function as Fourier transform pairs: 



1 

2^ 



P(k) k' 



sin ks 

ks 



d/c. 



(11) 



And finally, we compare these determinations with those 
obtained from the computation of the two point correlation 
function for each subsamplc using Equation 9. In Figure 4 we 
display the power spectra obtained for each subsamplc. The 
filled circles in this figure correspond to the power spectrum 
for the full sample of groups. As expected, the amplitude 
of the power spectrum increeises with the average mass of 
groups in the subsamples. If wc arc to estimate ^(s) using the 
determinations of P{k) presented in this section, it should 
be taken into account that a reliable computation requires 
knowledge of the power spectrum over a wide range of k 
values. Since we only obtain P{k) for k < 0.4 h Mpc~^, 
wo extend our estimates to larger wavcnumbcrs assuming a 
power law behaviour P{k) oc k~^ . Consequently, we fit this 
power law to each subsample using values of P{k) in the 
range 0.16 ^ A; ^ 0.4 h Mpc~^. The resulting power laws 
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With these extensions to large wavenumbers the two point 
correlation functions obtained applying Equation 11 to the 
results from each sample are plotted with solid lines in Fig- 
ure 5. In this figure, the filled circles correspond to the direct 
estimate of ^(s) for each subsample using Equation 9. The 
error bars are computed using the set of 10 2dFGGC mock 
catalogues described in section 3. From this plot we observe 
that there is a good agreement between both methods of 
estimating £,{s). The dashed lines in each panel correspond 
to the best fit obtained for ^(s) assuming a power law shape 
(eq. 10). The best fitting parameters obtained for each sub- 
samplc arc shown in table 1. 

Using the estimates of sq for each subsamplc, we pro- 
ceed to study the so — dc relation. The mean intergroup sepa- 
rations, dc, are computed using the analytical mass function 
prediction based on the ellipsoidal collapse model of overden- 
sitics derived by Slieth & Tormcn (1999) which is a very good 
fit to the 2dFGGC mass function measured by Martinez et 
al. (2002b). We compute the abundance n(> M.) of these 
systems in the ranges of masses defined for each subsample 
and estimate dc = n~^^^ . 

The resulting values of mean inter-group separation for 
the 2dFGGC subsamples are shown in table 1. The corre- 
sponding (dc. So) pairs are shown in Figure 6 by filled circles. 
For comparison, wc also plot previous determinations from 
groups in the Updated Zwicky Catalogue (GUZC) derived 
by Merchan, Maia & Lambas (2000) which shown a very 
similar behaviour with our results. The so — dc pairs shown 
in the range of scales corresponding to clusters of galaxies 
arc those obtained from the APM cluster survey by Croft 
et al. (1997, open squares). In previous works, some authors 
(Bahcall & West 1992, Bahcall & Cen 1992) have argued 
that Abell clusters are consistent with a linear so — dc re- 
lation described by so = 0.4 dc. This scaling law is plotted 
as a dashed line in Figure 6 and does not seem to provide a 
good description either of cluster data or of groups of galax- 
ies. This result was also reached by Croft et al. (1997) where 
they also computed the so — dc relation for clusters of galax- 
ies in a ACDM N-body simulation. Their results are repre- 
sented by crosses in Figure 6. As can be seen, our estimate of 
the So — dc relation for groups is an extension of the relation 
obtained by Croft et al. (1997) from numerical simulations 
from a ACDM model. Similar trends in sq — dc were also 
found in other CDM models, but with slightly different am- 
plitudes of correlation length (see Governato et al. 1999). 
We find that a simple fit of the form 



So = 4.7 dc 



(12) 



provides a satisfactory empirical description of the observa- 
tional data as well as the N-body simulation results. This 
empirical scaling law is shown by the solid line in Figure 6. 
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5 REDSHIFT SPACE DISTORTIONS 

In this section we study the redshift-space correlation func- 
tion calculated in the directions parallel and perpendicular 
to the line of sight, ^{a,n). This approach has led to the 
quantification of characteristics of the redshift-space distri- 
bution of galaxies and clusters of galaxies, such as the "fin- 
gers of God" , which are elongated structures seen in redshift 
surveys, originating from the random motions of galaxies 
inside clusters. The statistics that can be derived from the 
anisotropics in ^(cr, tt) are presented in section 5.1 and the 
results found for the diff'crcnt 2dFGGC subsamplcs can be 
found in section 5.2. Section 5.3 contains the results for pair- 
wise velocities corresponding to the different subsamples, 
and 5.4 the corresponding bias factors. 

5.1 Statistics from the anisotropies in the 
redshift-space correlation ftinction 

The quantification of the "fingers of God" effect comes 
from the measurement of the pairwisc velocities of galax- 
ies, which can be done by comparing a theoretical expres- 
sion for the correlation function, ^^{a, it), with the measured 
values, ^°(cr, tt) (the index o indicates quantities measured 
from observational data, and the index p denotes predicted 
quantities). Specifically, we compare curves of equal corre- 
lation function amplitude {^■^^'^) parametrised by r'^fix{9) 
and r^f,^ (d), where 9 is the angle measured from the direc- 
tion perpendicular to the lino of sight (Padilla et al. 2001), 
such that Ci'^r,'n^r) - S,^'''^ = 0, where = r'^ji^^id) cos{9), 
■Wr = rcfii,{9)sm{9), and the index x indicates either ob- 
served (o) or predicted quantities (p). 

The predicted correlation function is obtained by con- 
volving the real-space correlation function, ^(r) with the 
pairwisc velocity distribution function, f{w), following Bean 
et al. (1983). In order to do this, we calculate 

/DO 
[l + arMw -w'{r,r')]dw', (13) 
-oo 

where = r'^ + cr'^, and Ho is the Hubble constant, r' = 
n — w' /Ho (the prime denotes the line-of-sight component 

of a vector quantity) and 

w'{r,r') ~ -HoK{r){l + £,{r))-\-', (14) 

is the mean streaming velocity of galaxies at separation r. 

Following the usual prescription, we calculate the best- 
fit rms peculiar velocity, {w^)^^^ , for an exponential distri- 
bution. 

We adopted this pairwise velocity distribution as it has 
shown to be the most accurate fit to the results from nu- 
merical simulations (Ratcliffe et al. 1998, Padilla & Baugh 
2002). 

In this work we use a non-linear CDM correlation func- 
tion, for which we assume a given cosmology. In order to 
do this, we follow the recent results from satellite and bal- 
loon borne cosmic microwave background experiments such 
as WMAP (Spergel et al. 2003) and Boomerang (Netterfield 
et al. 2002) , and assume a ACDM cosmology with mass den- 
sity f2o = 0.28, baryon density Q.b = 0.05, vacuum energy 



density JIa = 0.67, and a CDM shape parameter F = 0.18, 
and (78 = 0.9. These parameters are compatible with con- 
straints from CMB data. Even though we include explicitly 
a baryon fraction in our models, the high level of noise of our 
observational data allow us only to measure the damping in 
power produced by baryons. This means that our choice of 
r implicitly allows for a combination of baryon fraction and 
a real value of F as discussed in Eisenstein & Hu (1996). We 
obtain the real-space correlation function by Fourier trans- 
forming the CDM power spectrum 

= ^^^/°°nfc)^fe^dfe, (16) 

and then use this to evaluate the theoretical prediction for 
^''(cr, tt) using equation (13). Here we assume a constant, 
scale independent, bias between the distribution of groups 
and mass, which has been shown to be a good approximation 
for groups and clusters of galaxies (Padilla & Baugh, 2001). 
Furthermore, by looking at figure 3, it can be seen that 
the values of correlation function slopes, 7, for the different 
samples are not significantly different (less than I.Sct away); 
also, figure 2 shows that a constant relative bias is a good 
approximation up to scales k < 0.45/i/Mpc, equivalent to 
X ~ lQh~ IMpc, where our analyses take place. 

We search for the optimum values of the scale inde- 
pendent bias parameter h and (w^)^''^ by minimising the 
quantity x'({w')"'', &), 

x\{wY^\ b) = Yyis.. m - r^,., {9i)f, (17) 

i 

where wo have chosen to compare a set of discrete levels, 
(^l^^ = 0.6, 0.8, 1.0, 1.2 and 1.4, of the redshift-space correla- 
tion function ^(cr, tt) amplitude (Padilla et al. 2001) instead 
of comparing values of the correlation function on a grid of 
a and tt distances as done previously in other works (see for 
instance Ratcliffe et al. 1998). 

5.2 Measurements of 2dfGGC ^(tr, tt) 

Using the estimator represented by equation (9), we measure 
the redshift-space correlation function of the 2dFGGC for 
the four subsamples studied in this paper. In this case, the 
quantities DD^ RR, and DR depend upon the separations 
parallel and perpendicular to the line of sight. The results 
can be seen in Figure 7 where the different panels show the 
correlation function for subsamples 1,2,3 and 4 (left to right 
hand side panels respectively). 

Concentrating on the anisotropics we observe in this 
figure, we find that the infall pattern expected from gravi- 
tational instability is seen in most of the group subsamples 
with the exception of subsample 3. Here we should bear in 
mind the increasingly noisy results that are obtained for the 
samples of higher virial masses. This is a simple consequence 
of the ever smaller number of more massive objects that adds 
a Poisson error of increasing magnitude to the estimate of 
the correlation function. Still, subsamples 1, 2 and 4 show 
the infall pattern for many contour levels. 

By inspecting the correlation lengths depicted by the 
thick solid lines in the panels of Figure 7, the increase in cor- 
relation length for subsamples with higher low virial mass 
limits is noticeable. This is in agreement with the direct mea^- 
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o-/(h"^ Mpo) cr/(h"' Mpc) o/(tC^ Mpc) o-/(h"' Mpc) 



Figure 7. The 2-point correlation function of the four subsamples of 2dF groups studied in this paper, in the coordinates a and 
TT (subsamples 1 to 4 from left to right hand side panels). The transitions between different shadings correspond to fixed values of 
^ = 0.6,0.8, 1.0, 1.2, and 1.4 levels, which are the values used later to infer the pairwise velocities, (w^)^/^. The thick line corresponds to 
the ^ = 1 level, and the thin lines show the expected contours in the absence of peculiar velocities. 



Table 1. Statistical results for the 2dFGGC. The first columns show several characteristics of the different 2dFGGC samples, including 
the minimum virial mass of the groups in each sample, number of groups, their average distance, the average immber of galaxy members 
identified per group, and the average internal velocity dispersion as measured from the identified member galaxies. Columns 7 and 8 
show the correlation function fit parameters, column 9 contains the values of mean inter-group separation, and the last two columns, the 
corrected values of group pairwise velocity and bias factors obtained from the redshift-space distortions analysis. 



Sample 


•A4vir min. Ngrp {^^) '^member 


so 


7 




4 




[H-^Mq] [km s-l] [km s-i] 


[h-iMpc] 




[ft-iMpc] [km s-i] 





All 2198 265 ± 150 8.9 ± 0.3 -1.6 ±0.1 1.82 ± 0.36 

Sample 1 5 x lO^^ 1996 30333 6.9 ± 4.6 283 ± 149 9.6 ± 0.4 -1.7 ±0.1 8.94 ZSOI^Jq 1.92 ± 0.38 

Sample 2 2 x lO^^ 1448 33035 7.5 ± 5.0 333 ± 143 10.8 ± 0.6 -1.7 ±0.2 15.00 395+95 2.04 ± 0.41 

Sample 3 5 x lO" 920 36336 8.2 ± 5.8 396 ± 141 11.7 ± 0.6 -2.0 ± 0.2 20.93 540l*oo 2.24 ± 0.45 

Sample 4 1 x lO" 540 38631 8.6 ± 6.4 464 ± 142 13.4 ± 0.9 -2.0 ± 0.2 28.93 49515^5 ^.51 ± 0.50 



surements of redshift-space correlation function performed 
in section 4.2. 



5.3 Pairwise velocities and correlation lengths in 
observationEil cluster samples 

In order to summarise the results obtained from the study 
of the anisotropies in the correlation function of 2dFGGC 
subsamples, we present the values of pairwise velocities ob- 
tained from the correlation functions shown in Figure 7. 

We find the group pairwise velocities by minimising 
equation 17. The real space correlation function used in this 
equation corresponds to a CDM power spectrum with spec- 
tral index Ua = 1.0, and parameters flo = 0.28, Oa = 0.67, 
a baryon density Ob = 0.05, a CDM shape parameter 
r = 0.18, and as = 0.9. 

The results for the pairwise pairwise velocities found 
from the correlation functions of the four 2dFGGC subsam- 
ples studied here axe presented in Figure 8. The errorbars 
correspond to a 1 — cr standard deviation, calculated using 
5 estimates of pairwise velocities obtained from 5 different 
correlation function level contours. Wo note that this is at 
best a lower limit to the uncertainty in (w^)^^^, since the 
results from different ^(fT, tt) can be highly correlated. The 
ordering of samples is such that an increase in the x direc- 



tion corresponds to an increase in the sample's lower virial 
mass limit. 

The assumption that structures form by hierarchical 
clustering implies that the values of pairwise velocities found 
for these subsamples should have been smaller for higher 
mass limits, since more massive objects take longer to form 
virialized structures. This trend is expected in the direct 
measurements of pairwise velocities of groups identified from 
a ACDM simulation. In order to test this hypothesis we take 
20 subsets of groups with the same mass limits and number 
of members as the 2dFGGG subsamples from the simula- 
tion, and measure their 1-dimensional pairwise velocities. 
The downward hatched region in Figure 8 corresponds to 
the 1-cr confidence levels for the group pairwise velocities 
from the numerical simulations. Even though the trend of 
decreasing pairwise velocities for higher masses is visible in 
the simulation results, the errors make it also consistent with 
a constant value. As can bo soon, the agreement between 
subsample 1 and the simulations is remarkable. liowcver, 
more massive groups in subsamples 2, 3 and 4 show higher 
values of (w^)^''^. This behaviour is in the opposite direction 
to what we expected from hierarchical clustering arguments. 
A possible cause for this disagreement is the likely inclusion 
of spurious groups affecting the low mass sample. However, 
as shown by Padilla & Lambas (2003a), this would produce 
an enhancement of the pairwise velocities, and therefore our 
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and (T„ in Table 1). An expression for this uncertainty as 
a consequence of the finite number of galaxy distance mear- 
surements per group is, 



A{cz) = 



2 3 
Sample Number 



Figure 8. Pairwise velocities obtained for the different 2dFGGC 
subsamples (open symbols with errorbars). The samples are dis- 
tributed along the x-axis, with samples of increasing lower virial 
mass limit towards the right. Errorbars are obtained from the es- 
timates of pairwise velocities at 5 different correlation function 
contour levels. The solid line shows the corrected pairwise ve- 
locities obtained by subtracting in quadrature the error in the 
distance measurement of groups, and the corresponding 1-<t con- 
fidence level is represented by the upward hatched region. The 
downward hatched region shows the I-ct confidence level obtained 
from 20 subsamples of groups identified from a ACDM simulation. 
Each of these subsamples contain the same number of groups as 
the 2dFGGC subsamples. 



results for sample 1 would be overestimated, making the dis- 
agreement more severe. However, Padilla & Lambas (2003a) 
also show that in the case of groups identified from red- 
shift surveys, this effect is not important. Another possi- 
ble factor influencing our results is the increasing values of 
the average distance to groups in more massive subsamples, 
which is brought by selection function effects. The average 
group distance in subsamples 1,2,3, and 4, is larger for the 
most massive samples, as can be seen in Table 1. It is pos- 
sible that the increase in average distance translates into a 
larger distance measurement error, which in turn tends to 
enhance the anisotropies in the correlation function along 
the line of sight, thus explaining the results obtained here. It 
should be noted though, that the amplitude of errors needed 
to obtain agreement with the simulations is small. For in- 
stance, for samples 3 and 4, a distance measurement error 
of A{w^)^^^ ~ 450km/s is enough to explain the discrepan- 
cies between the results from observations and the numerical 
simulation. This error is equivalent to a 1% uncertainty in 
the group distance measurements at the mean depth of these 
subsamples (see Table 1). 

The origin of the uncertainty in the group distance mea- 
surement can be related to the fact that the number of 
galaxies {rimember) uscd in the distance determination of the 
groups shows little change among the subsamples, whereas 
the internal velocity dispersion of galaxies in groups (cr„) 
shows a clear enhancement (see average values of rimember 



rimember 



(18) 



wliieli is based on the assumption that these galaxies were 
randomly selected from a gaussian distribution in redshift- 
space, of width (t„. Using the average rimember in the group 
subsamples, we subtract yJ\2)A.{cz) in quadrature from the 

estimate of (w^)^^^, and show the corrected {nP)V^ in a 
solid line in Figure 8 as a function of subsample number. 
The upward hatched region shows the range covered by a 1- 
(T uncertainty in {w^)^^, as it results from considering the 
widths of the distributions in rimember and a„ for the in- 
dividual group subsamples, and the original uncertainty in 
(w^)^''^ (the widths in rimember and cr^ are given in Table 
1 as the uncertainty in these quantities). As can be seen, 
the corrected values are only slightly smaller than the orig- 
inal {vP)^^'^. However, the amplitude of the uncertainties 
is noticeably increased, improving in most of the cases the 
agreement with the ACDM simulation. This is not enough 
to obtain an agreement for subsample 3, though, which had 
also failed in showing the expected infall pattern in figure 7. 

We remark however that the assumed value of erg also 
influences our results on {w^)^^^ , but to a lesser extent. As- 
suming (Ts = 0.7 lowers our estimates of pairwise velocity 
by ~ 10%, but leaves the uncertainty unchanged. However, 
we bear in mind the fact that these uncertainties can be 
underestimated due to correlations between the individual 
estimates of {w^)^^^, as explained above. 

Another important point to be noticed is that even for 
sample 3 (our sample with no infall pattern), the obtained 
value of corrected group pairwise velocity (w^)^^^ remains 
one of the lowest values found for systems of galaxies (see 
Padilla & Lambas 2003b for a comprehensive analysis of 
several cluster samples), and is still in agreement with the 
pairwise velocities found for galaxy catalogues {w'^)^^^ ~ 
(450 ± 100)km/s (Ratcliffe et al. 1998, Padilla et al. 2001). 



5.4 Bias factors 

The comparison between the measured and predicted 
redshift-space correlation functions, as explained in sec- 
tion 2, involves the real-space bias (6^) between the group 
and mass distributions, where the latter is described by a 
CDM model. In this work, we measure by minimising 
x'^{{w^y^'^,bi){eq. 17), where (w^)^^"^ corresponds to the 
group pairwise velocities found in the previous section. As 
mentioned above, depends upon the real space correla- 
tion function which is obtained from a CDM power spectrum 
with a fixed set of cosmological parameters (see previous sec- 
tion). We minimise x^i and show the resulting bias factors 
in Figure 9 (open circles with errorbars) where the x-axis 
shows the 2dFGGC subsamples in order of increasing group 
mass. The errors on the bias parameter come from the re- 
sults of minimising equation (17) for 5 different levels of the 
correlation function. We also include a 20% error added in 
quadrature to allow for an uncertainty in the assumed value 
of as = 0.9 (for instance, assuming as = 0.7 would produce 
an enhancement of a 13% in our estimates of bias factors). 
The values of bias factors and their uncertainties can be 
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found in Table 1, where we also show the result obtained for 
the full 2dFGGC sample. 

These values can be compared to what is expected from 
a particular CDM model, using the expression for an ef- 
fective CDM bias as presented in Padilla & Baugh (2002), 
which is a weighted average of the bias of a sample of clus- 
ters of mass M, derived by Sheth, Mo & Tormen (2001), 
and requires knowledge of the space density of the cluster 
sample. The results from using this equation are valid if the 
cluster sample is complete above a minimum mass threshold. 
The main source of systematics arising from this assumption 
comes from a possible underestimation of the space density 
induced by the incompleteness of the group sample. How- 
ever, a lower limit for the completeness of the 2dFGGC sub- 
samples is ~ 0.7, and would only induce a change of ~ 5% 
in the CDM bieis factor. 

The hatched region in figure 9 shows the range of ex- 
pected values of a CDM bias factor for the same space den- 
sities present in the 2dFGGC subsamples, assuming 0.75 < 
as < 0.9 (a larger as corresponds to a lower CDM bias fac- 
tor). As can be seen, the 2dFGGC groups bias factor tends to 
increase with the minimum virial mass of the sample as ex- 
pected from hierarchical clustering. Furthermore, the CDM 
bias is within the 1 — a confidence levels around the mear 
sured bias factors for most of the samples, the exemption 
being the least massive subsample, for which the expected 
CDM bias factor is smaller. We investigated the effects of 
an underestimation of the group pairwise velocities on the 
value of the bias, and found that an increase of a factor ~ 2 
in (ui^)^^^ is enough for obtaining a result compatible with 
the CDM expectations. 

As mentioned above, there is a possible ~ 5% system- 
atic error in the CDM bias factors. In order to visualise the 
effects of such an error, we show two solid lines in figure 9 
which limit the range of CDM bias factors after correcting 
for this systematic error. As can be seen, this does not sig- 
nificantly affect the comparison between the 2dFGGC and 
CDM model biases. 

The bias factors calculated in this section correspond 
to real space quantities and are relative to the distribu- 
tion of mass which has an assumed amplitude of density 
fluctuations erg = 0.9. On the other hand, the relative 
bias found from the ratio between the galaxy and group 
power spectra in section 4.1, bs{k), corresponds to a red- 
shift space quantity, and is relative to the distribution of 
galaxies. The fact that this value corresponds to redshift- 
space, allows us to make a rough estimate of what would 
the corresponding real-space value be, based on the studies 
of Padilla & Baugh (2002), who find that the ratio between 
real and redshift space biases for cluster samples is roughly 
br/bs = 1.25. Assuming that this relation can be also ap- 
plied to the 2dFGGC, the redshift-space bias factor found 
from the power spectrum analysis corresponds to a real- 
space br{k) = 1.9 ±0.4, in the range of scales studied in this 
paper, 0.025 < k/fiMpc^^ < 0.45, consistent with our esti- 
mates of the real-space bias parameter from the anysotropies 
in the correlation function. Assuming a different value of 
as = 0.7 produces no significant changes in this result. 



: T 

2.5 - ) 



1.5 - 




1 2 3 

Sample Number 



Figure 9. Bias factors for the 2dFGGC subsamples studied in 
this paper. Tlie circles with errorbars show the results from min- 
imising equation 17, and the hatched region shows the range of 
predicted effective bias from CDM and the Sheth, Mo & Tormen 
(2001) mass function for 0.75 < as < 0.9. The solid lines limit 
the same range, after correcting for a systematic error ~ 5%. The 
different samples are displayed along the x-axis, ordered in such 
a way as to have higher virial masses towards the right. 



6 CONCLUSIONS 

We analysed the 2dFGGC constructed by Merchan & Zandi- 
varez (2002), one of the largest group catalogues constructed 
at present. Our analysis considers galaxy groups as tracers 
of the large scale structure in the universe. The statisti- 
cal tools adopted for this work are the power spectrum of 
density fluctuations and the two point correlation function, 
both calculated using redshift-space data. We calculated the 
group power spectrum and found that its shape is very simi- 
lar to that of the galaxy power spectrum, but shows a higher 
amplitude as expected for higher mass systems in hierarchi- 
cal clustering. The measurement of a relative bias between 
groups and galaxies in redshift space quantifies this differ- 
ence in amplitude, and results in an almost constant relative 
bias value of bs{k) ~ 1.5 on the scales probed by our anal- 
ysis, 0.025 <k/h Mpc"^ < 0.45. 

The estimate of the two point correlation function for 
the group sample is found to be well fitted by a power law 
of the form ^ — (s/so)^ with a correlation length so = 8.9± 
0.3 /i"^Mpc and a slope 7 = -1.6 ± 0.1. 

We also analyse the variation of the 2dFGGC groups 
clustering with group mass by studying the so — dc relation. 
For this purpose, we split the group sample in four sub- 
samples according to difi'erent ranges of group virial masses. 
Firstly, we check that the computation of the correlation 
function for each subsample using the direct method defined 
by Landy & Szalay (1993), yields a remarkably similar result 
to that obtained by Fourier transforming the corresponding 
power spectrum estimate. Secondly, the obtained correlation 
length So tends to increase when using more massive samples 
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as previously observed by Merchan, Maia & Lambas (2000, 
see Table 1). 

From the computation of the mean intergroup distance 
dc, the resulting so — dc relation is quite similar to that 
obtained by Merchan, Maia & Lambas (1993) for the GUZC. 
If we extend this relation with the so — dc relation obtained 
for cluster of galaxies by Croft et al. (1997) we observe that 
the universal scaling law of Bahcall & West (1992), so = 
0.4 dc, provides a poor description to the data. We find 
that a fit described hy the law sq = 4.7 dc""^^, is capable 
of describing both, the observational so — dc relation found 
for groups and clusters, as well as the results from N-body 
numerical simulations of a A CDM model presented by Croft 
et al. (1997, see also Govcrnato et al. 1999). 

The results from the study of the anisotropics in the 
redshift space correlation function show that in general, the 
groups in the 2dFGGC present the expected infall pattern 
indicating that these systems have not yet formed virialized 
structures themselves. This can be seen from the correlation 
function contours of sample 1, which includes the majority 
of the 2dFGGC groups. 

The study of the redshift-space correlation functions 
also allowed us to make an estimate of the groups bias factor 
relative to the distribution of the mass. In order to achieve 
this, it was necessary to assume a CDM model with a fixed 
value of (78 = 0.9. We compare our estimates of bias factors 
to what is expected for such a cosmological model, where the 
mass function is well described by the Shotli, Mo & Tormcn 
(2001) model. We find very good agreement for the samples 
2, 3 and 4. 

We calculate the group pairwiso velocities for the dif- 
ferent 2dFGGC subsamplcs, and correct them from group 
distance uncertainty effects ((w^)^^, see table 1). We note 
that the pairwise velocity measured for samples 1 and 2, 
= (280l5;o)km/s and (w^)^^ = (395l^^)km/s, are 
noticeably smaller than that found for galaxies, {tu^)^^^ = 
(450 ± 100)km/s, in agreement with gravitational instabil- 
ity expectations. This agreement comes not only from a 
group pairwise velocity which is qualitatively smaller than 
the measurements found for galaxies, but also from a com- 
parison with measurements from actual simulated groups, 
embedded in one of the most readily motivated cosmologi- 
cal scenarios, a ACDM model. 
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